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ABSTRACT 

High-resolution spectra of the hot white dwarf G191-B2B, covering the wavelength region 905-1187A, 
were obtained with the Far Ultraviolet Spectroscopic Explorer (FUSE). This data was used in conjunction 
with existing high-resolution Hubble Space Telescope STIS observations to evaluate the total H i, D i, O i 
and N i column densities along the line of sight. Previous determinations of N{D i) based upon GHRS 
and STIS observations were controversial due to the saturated strength of the D i Lyman a line. In the 
present analysis the column density of D i has been measured using only the unsaturated Lyman /3 and 
Lyman 7 lines observed by FUSE. A careful inspection of possible systematic uncertainties tied to the 
modeling of the stellar continuum or to the uncertainties in the FUSE instrumental characteristics has 
been performed. The column densities derived are: log V(D i) = 13.40 ±0.07, log V(0 i) = 14.86 ±0.07, 
and log V(N i) = 13.87 ± 0.07 quoted with 2cr uncertainties. 

The measurement of the H i column density by profile fitting of the Lyman a line has been found to be 
unsecure. If additional weak hot interstellar components are added to the three detected clouds along the 
line of sight, the H I column density can be reduced quite significantly, even though the signal-to-noise 
ratio and spectral resolution at Lyman a are excellent. The new estimate of N(H i) toward G191-B2B 
reads: log V(H i) = 18.18 ± 0.18 (2cr uncertainty), so that the average (D/H) ratio on the line of sight 
is: (D/H)= l.GG^o'e ^ uncertainty). 

Subject headings: ISM: abundances : ISM — ultraviolet 


1. INTRODUCTION 

Deuterium is thought to be produced in significant 
amount only during primordial Big Bang nucleosynthesis 
(BBN), and to be thoroughly destroyed in stellar interi¬ 
ors. Deuterium is thus a key element in cosmology and in 
Galactic chemical evolution (see e.g. Audouze & Tinsley 
1976; Gautier & Owen 1983; Vidal-Madjar & Gry 1984; 
Boesgaard & Steigman 1985; Olive et al. 1990; Pagel 1992; 
Vangioni-Flam & Casse 1994; Vangioni-Flam et al. 1995; 
Prantzos 1996; Scully et al. 1997; Casse & Vangioni-Flam 
1998; Tosi et al. 1998). Indeed, its primordial abundance 


is the best tracer of the baryonic density parameter of the 
Universe, fls, and the decrease in its abundance during 
galactic evolution traces the cosmic star formation rate at 
various epochs. 

The first (indirect) measurement of the deuterium abun¬ 
dance was carried out using ^He in the solar wind, giv¬ 
ing the presolar value D/H~ 2.5 ± 1.0 x 10“® (Geiss & 
Reeves 1972). The first measurements of the D/H ratio 
in the interstellar medium (ISM) were reported shortly 
thereafter by Rogerson & York (1973), and their value 
(D/H)~ 1.4±0.2 X 10“^ has remained a landmark average 
value for the interstellar D/H ratio. Finally direct mea- 
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surements of the primordial (D/H) ratio in low-metallicity 
material at high redshift have been successfully carried 
out these past few years (e.g., Buries 2001 for a review, 
and references therein). The values derived cluster around 
(D/H)^ 3 X 10 ^ although with significant dispersion, 
which may or may not be real. Quite similarly the mea¬ 
surements of the (D/H) ratio in the Galactic ISM towards 
hot stars with the Copernicus satellite lead to many eval¬ 
uations of D/H (see e.g. York and Rogerson 1976; Vidal- 
Madjar et al. 1977; Laurent et at 1979; Ferlet et al. 
1980; York 1983; Allen et al. 1992) which also show dis¬ 
persion around the above York & Rogerson (1973) value. 
This dispersion has been recently confirmed by IMAPS 
observations (Jenkins et al. 1999; Sonneborn et al. 2000), 
indicating that the D/H ratio may vary by a factor ~ 3 in 
the solar neighborhood, i.e., within a few hundred parsecs. 

In this paper we present a new determination of the D/H 
ratio on the line of sight to the nearby DA white dwarf 
(WD) G191-B2B based on observations obtained with the 
Far Ultraviolet Spectroscopic Explorer {FUSE, Moos et al. 
2000; Sahnow et al. 2000). This paper is one of a series in 
this volume describing the first results of the FUSE (D/H) 
program in the Local ISM (LISM). This program and its 
results are summarized in the overview paper by Moos et 
al. ( 2002 ). 

Observing white dwarfs has many advantages over hot 
and cool stars for studying the D/H ratio, as explained in 
Vidal-Madjar et al. (1998): these targets can be chosen 
close to the Sun, in order to avoid a complex line of sight 
structure, and in the high temperature range, so that the 
interstellar absorption is superimposed on a smooth stellar 
continuum. The risk of contamination by low column den¬ 
sity H I fluffs possibly present in the hot star winds (Gry, 
Lamers & Vidal-Madjar 1984) is negligible for WDs, and 
their hot continuum offers the possibility of observing the 
numerous UV lines of N i and especially O i, which is a 
very useful tracer of H i and D i. The (D/H) ratio has al¬ 
ready been measured toward four white dwarfs, using the 
HST: G191-B2B (Lemoine et al. 1996; Vidal-Madjar et al. 
1998; Sahu et al. 1999), HZ43A (Landsman et al. 1996), 
Sirius B (Hebrard et al. 1999) and Feige 24 (Vennes et al. 
2000). For HZ43A, Sirius B and Feige 24 the average D/H 
values obtained are compatible with the local ISM (LISM) 
D/H determination (Linsky 1998) made within the Local 
Interstellar Gloud (the LIG in which the sun is embed¬ 
ded), although in the case of Sirius B this compatibility is 
marginal. 

In the case of GI9I-B2B, it was found that the line 
of sight comprises one neutral region corresponding to 
the LIG, and two more ionized absorbing components 
(Lemoine et al. 1996; Vidal-Madjar et al. 1998). The 
average (D/H) ratio (defined as the ratio of the total col¬ 
umn densities of D i and H i) was found to be (D/H)= 
1.12 ± 0.08 X 10“^ (Vidal-Madjar et al. 1998), signifi¬ 
cantly different from the value measured toward Gapella 
(D/H)lic = 1.5 ±0.1 X 10-5 (Linsky 1998). The (D/H) 
ratio measured toward G191-B2B has been contested by 
Sahu et al. (1999), who used STIS data and concluded to 
the presence of two interstellar components only, and to a 
D/H ratio compatible with that observed toward Gapella. 
The disagreement resides in the evaluation of the total 
D I column density (Vidal-Madjar 2000; Sahu 2000), and 


arises presumably because the Lyman a D line is satu¬ 
rated and the column density is thus sensitive to the line 
profile. The number of components assumed on the line 
of sight may also introduce differences between the analy¬ 
ses of these groups (see Vidal-Madjar 2001 for a detailed 
discussion). 

In the present work we re-examine these issues, making 
use of high quality FUSE and STIS observations of G191- 
B2B. We first measure the N(D i), N{0 i) and N(N i) col¬ 
umn densities using the unsaturated lines of these elements 
in the FUSE datasets, notably Lyman (3 and Lyman 7 for 
D I (Section 2). We then analyze the recent high qual¬ 
ity STIS observations in Section 3, and provide explicit 
evidence for the presence of three absorbing components 
(at least) on the line of sight. We also provide a refined 
estimate of the total H i column density. All throughout 
this work, considerable effort has been put on quantifying 
possible systematic uncertainties related to fixed-pattern 
noise, detector artifacts, background uncertainties, wave¬ 
length calibration and modeling of the stellar continua, 
as well as to the velocity structure of the line of sight. 
In particular, we argue in Section 3.3 that previous esti¬ 
mates of the total V(H i) are subject to a large systematic 
uncertainty related to the possible presence of additional 
weak [V(H i)< 10^^ cm”^] hot {T ~ lO^K) components. 
This effect may have a large impact on our understanding 
of the observed variations of the (D/H) ratio in the ISM, 
as it may affect other lines of sight, and is the subject 
of a companion paper (Vidal-Madjar & Ferlet 2002). We 
provide a summary of our results and a short discussion 
in Section 4; an overall discussion of the FUSE results is 
given by Moos et al. (2002). 

2. COLUMN DENSITIES OF D I, O I AND N I MEASURED 
WITH fuse 

The FUSE instrument (Moos et al. 2000; Sahnow et 
al. 2000) gives access to the wavelength range 905-II87A 
on eight detector segments, through three different en¬ 
trance apertures (the large, medium and narrow apertures 
- LWRS, MDRS and HIRS, respectively) and four chan¬ 
nels (two LiF and two SiC). In particular, it gives access to 
the higher order unsaturated Lyman lines of D i and can 
thus offer a reliable estimate of N(D i). Similarly, many 
weak lines of O i and N i are present in this bandpass, 
as shown in Table 1. The FUSE spectrum of G191-B2B 
in the SiClB channel observed with the HIRS aperture is 
shown in Fig. 1, in which the lines of D i, N i, O i used in 
our study are indicated. 

The white dwarf G191-B2B has been observed several 
times with FUSE through the three different entrance 
apertures; the observation log is given in Table 2. Since 
the target is point-like, the resolutions obtained with these 
apertures should not be too different, but there is a clear 
variation due to some off-axis effect when moving from 
one aperture to another. The main problem comes from 
the geocoronal Lyman (3 emission whose strength increases 
with the surface area of the aperture, and which blends 
with the interstellar H i absorption. This leads us to em¬ 
ploy different approaches in the analysis of the data, de¬ 
pending upon the slit used, as discussed below. 

The observations were treated through the version 1.8.7 
of the CALFUSE pipeline. The data were then collected 



3 


channel by channel by series of subexposures. Due to the 
sensitivity of FUSE and the brightness of the source, the 
data were obtained in the HISTOGRAM mode (a spectral 
imaging mode, in which a two-dimensional spectrum was 
accumulated onboard). Each subexposure has a S/N ratio 
that is high enough to clearly see the strong interstellar or 
photospheric features. The different subexposures could 
then be easily aligned on top of each other to compensate 
for the slow thermal drifts that displace the wavelength 
scale from one subexposure to another. These drifts never 
translates into more than ±5 pixels. A high S/N spectrum 
is thus reconstructed in each channel and for each aper¬ 
ture. This process preserves the ultimate FUSE resolution 
and it also partly eliminates some of the detector’s fixed 
pattern noise by acting like a random FP-SPLIT proce¬ 
dure. 

Since the EUSE sampling is ~ 10 pixels per resolution 
element, we rebinned the spectra by three. We measure 
column densities by profile fitting of the various lines ob¬ 
served, using the code Owens. f developed by one of us 
(M. L.) and the French EUSE team. Our general proce¬ 
dure for data analysis is as follows. We split each spectrum 
into a series of small sub-spectra centered on absorption 
lines to be analyzed, whose width is a few A depending 
on the line density and fit all lines contained in all sub¬ 
spectra simultaneously. Between different sub-spectra, the 
zero-point velocity offset is left to vary to compensate for 
the wavelength calibration ripple effects in FUSE data; the 
Line Spread Function is generally taken to be a Gaussian 
whose FWHM varies freely during the fitting procedure 
and independently from sub-spectrum to sub-spectrum. 
It is also checked that considering a double Gaussian LSF 
with free widths, free amplitude ratio and zero pixel offset 
between the two Gaussian components, does not change 
the results. The background levels used for each line are 
those evaluated at the bottom of the closest Lyman line 
(see Fig. 2 for an example), which leads to an upper limit 
for the derived column densities; fixing the background 
level at the 0 value gives access to the lower limit of the 
evaluated column densities. This uncertainty is taken into 
account in our procedure; however EUSE is operating in 
first order grating mounts and consequently the stray light 
level is generally quite low, as can be seen in Fig. 2 for 
Lyman 7 . This analysis procedure allows to reduce sys¬ 
tematic uncertainties tied to fixed pattern noise effects or 
detector artifacts and uncertain calibration, since different 
lines of the same elements observed through different aper¬ 
tures and in different channels are fitted simultaneously. 
For instance, if a given line is subject to an artifact in the 
dataset, the of the fit in the region around this line, 
or equivalently the quality of the fit of this particular line 
should stand out in the overall fit if the number of lines 
analyzed is sufficiently large, as is the case here. More¬ 
over, by letting the unknown instrumental parameters vary 
freely, i.e., the LSF shape and width, the background flux, 
the wavelength zero-point and the shape of the continuum, 
and be optimized simultaneously with the physical param¬ 
eters that define the absorbers, one marginalizes the final 
result over these unknowns. Therefore the errors attached 
to the uncertainty on these above characteristics are in¬ 
cluded in the final errors quoted for the column densities. 
Finally, all lines are also analyzed one by one, i.e. each 


independently of all others, in order to check for overall 
consistency. A detailed discussion of this procedure and 
survey of possible systematic errors tied to FUSE data is 
given in a companion paper by Hebrard et al. (2002). 

2.1. O i and N 1 

Since FUSE cannot resolve the velocity structure of the 
line of sight (see Section 3), we have used only unsaturated 
lines of D I, N I and O i and assumed only one interstellar 
component to be present. Therefore, all measured col¬ 
umn densities in this Section are integrated over the line 
of sight. It has been verified that using three interstel¬ 
lar components in the fit of the EUSE did not change the 
conclusions as long as only non-saturated lines were con¬ 
sidered; an example of single vs multiple component fits is 
shown in Fig. 2. 

Note that HST observations cannot give a reliable mea¬ 
surement of the O I column density without a complete 
knowledge of the velocity and broadening structure of the 
line of sight since the only absorption line available at 
I 3 O 2 A is strongly saturated. In contrast the N i triplet 
at 1200A in the GHRS or STIS bandpass is only moder¬ 
ately saturated and may provide a precise estimate of the 
total A^(N i). The various measurements of 7V(N i) and 
N{0 i) are listed in Table 3. 

Since the contamination of interstellar lines by weak 
photospheric lines is a likely possibility in the EUSE spec¬ 
tral range due to the high density of atomic lines in this 
bandpass, we excluded from the fit all lines that presented 
an obviously excessive equivalent width or column density 
relative to the others. As an example as shown on Fig¬ 
ure 3 the N I line at 954.1042A is very probably blended 
with a photospheric line which shows up at nearly the 
same wavelength in the NLTE spectrum calculated with 
the TLUSTY/SYNSPEG codes (described in Section 3). 
Such blending becomes particularly important in the case 
of G191-B2B because the total O i or N i column densities 
are small, and their absorptions are relatively weak fea¬ 
tures which may be more easily blended with undetected 
photospheric ones. This is much less important for the 
other FUSE (D/H) targets (see the synthesis of FUSE 
results by Moos et al. 2002, and the individual line of 
sights studies by Friedman et al. 2002, Hebrard et al. 2002, 
Lehner et al. 2002, Sonneborn et al. 2002 and Wood et 
al. 2002) which present larger column densities, except for 
the white dwarf HZ43A (Kruk et al. 2002) which presents 
a pure H atmosphere (Barstow, Holberg & Koester 1995; 
Dupuis et al. 1998) and thus no possible blend with photo¬ 
spheric lines. However note that this line selection process 
may induce in the case of G191-B2B a slight underestima¬ 
tion of the total O i or N i column densities. 

The O I and N i lines used in this study are listed in 
Table 1. Some of them were rejected either because they 
are saturated (marked “strong” in Table 1), or because 
of a blend with a photospheric feature (marked “blend”), 
or because they are too weak and subject to systematic 
effects of the noise (marked “weak”). Some were kept in 
one channel but not in an other if a nearby detector defect 
was identified. For this reason, having the same spectrum 
observed in different independent channels was extremely 
helpful. 

By following the general method depicted above, the 
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final column density determinations are : 

log N{0 i)tot = 14.86 (±0.05) (±0.04) 

log iV(N ijtot = 13.87 (±0.05) (±0.05) 

where the errors are first statistical {2a) and second sys¬ 
tematic (95.5% C.I.). The statistical errors were evaluated 
through the technique, an example of which is given 
for the N{D i) determination to follow. However, one im¬ 
portant comment is of order at this point. The total 
for all spectral lines Htted simultaneously (H i, D i, O i 
and N i) in the case of the HIRS observations is equal to 
3211.44 for 1843 degrees of freedom (d.o.f.), even though 
the fit appears very satisfactory by eye. This value lies 
well above the upper 95.5% confidence level limit for sta¬ 
tistical fluctuation at 1966.71. We interpret this large 
value as a result of a systematic underestimate of the mag¬ 
nitude of the individual pixel errors by the instrument. In 
order to quantify this uncertainty, we have completed low 
order polynomial fits over relatively smooth and flat con¬ 
tinuum sections close to each of the spectral lines analyzed. 
These simple fits gave a of 688 for 500 degrees of free¬ 
dom, again above the 95.5% c.l. upper limit at 563.9. This 
indicates that indeed the individual pixel errors are under¬ 
estimated by instrument, probably due to the presence of 
fixed pattern noise in the data. The discrepancy in terms 
of ratio of measured to expected y^ is not as strong as for 
the best-fit y^ values, but this may be due to differences 
in the error levels for pixels located on the continuum and 
those located at the bottom of the saturated Lyman lines. 

In order to compensate for these effects, we have thus 
decided to rescale all y^ measured, which amounts to an 
overall rescaling of the individual pixel error estimates. 
In particular, we have chosen to divide all y^ by 1.84 = 
3211.44/1741.64, where 1741.64 corresponds to the lowest 
expected y^ at the 95.5% conhdence level for 1843 d.o.f. 
This rescaling rescales (increases) the error bars on the 
physical parameters that we derive since the Ay^ contour 
levels are also rescaled. In effect, in order to derive the 3a 
error bar around a single parameter minimum value, we 
search the values of this parameter that lead to an increase 
in (non-rescaled) y^ of 9 x 1.84 = 16.6 instead of 9. Our 
choice of rescaling is thus very conservative with respect 
to the Hnal error bars measured. All presented FUSE sta¬ 
tistical errors will be evaluated following this method. 

The systematic errors quoted above for O i and N i re¬ 
flect the range of values obtained for the best-fit solutions 
for different apertures, given in Table 3. Note that in this 
table, the errors quoted include the satistical error plus 
part of the above systematic errors, since it accounts for 
the differences measured between different observations, 
different channels and different spectral ranges where lines 
of the species concerned are seen, but observed through the 
same aperture (see the discussion in Hebrard et al. 2002). 
Note that the total O i or N i column density evaluations 
made with the LWRS±MDRS are in both cases lower than 
those made with the other apertures. We do not know the 
cause of this apparent systematic effect, but included it in 
the overall uncertainty on our result. 

The quadratic sum of statistical and systematic uncer¬ 
tainties then leads to: 

log N{0 i)tot = 14.86 (±0.07) 

log A(N i)tot = 13.87 (±0.07) 


2.2. Measurement of N(D i) 

In the FUSE bandpass, D i is clearly detected at Ly¬ 
man /3, weakly at Lyman 7 , and not at Lyman d. We 
list in Table 4 the different estimates of the total column 
density of D i that have been obtained from HST observa¬ 
tions of Lyman a, and which we measure using the present 
FUSE data. Here as well, we assumed that only one com¬ 
ponent is present on the line of sight since the D i lines 
are optically thin. Again we checked the impact of a more 
complex line of sight structure and found it to be negligi¬ 
ble, as expected. We carried out the same investigations 
of systematic effects as for O i and N i, using the approach 
detailed above and in Hebrard et al. (2002). 

There is however one essential difference between D i 
and O I or N I in terms of systematics. Namely the con¬ 
tinuum to the D I absorption is provided by the blue wing 
of the corresponding H i absorption, while the continuum 
around the O i and N i is very smooth. It is thus necessary 
to measure the possible systematic tied to the estimate of 
the stellar continuum on the measurement of 7V(D i). In 
order to do so, we have measured iV(D i) with and without 
beforehand correction of the data by a theoretical NLTE 
stellar continuum (described in Section 3). Some corre¬ 
sponding fits are shown in Fig. 2, and the best-Ht to the 
HIRS data is shown in Fig. 4. The effect of any assump¬ 
tion on the stellar continuum was found to be negligible 
as compared to the statistical uncertainty. 

However, we have found that the value of N(D i) mea¬ 
sured if the H I lines are excluded from the Ht and their 
blue wing modeled by a polynomial is higher than when 
the H I are included in the fit. This systematic effect 
may result from the following: in the latter (standard) ap¬ 
proach, where both D i and H i are fitted together, the 
position of the D i absorption is tied to that of H i, i.e., 
they have the same radial velocity and their zero point 
wavelength offsets are the same as both lines appear in the 
same sub-spectrum. However, in the former approach, the 
velocity of D i is not subject to this constraint. These two 
approaches are denoted in Table 4 as “D i & H i” (stan¬ 
dard approach) and “D i, no H i” for the other one. One 
should note that a discrepancy in radial velocity between 
the D I line and its counterpart H i could be attributed to 
small scale ripples in the FUSE wavelength calibration or 
to the presence of weak H i absorbers which would shift 
the position of the H i line with respect to that of the D i 
line in which these absorbers would not be felt. In this 
respect the approach in which D i is fitted independently 
of H I seems more adequate; however the interpolation of 
the blue wing of the H i line by a polynomial is not always 
straightforward. We should mention that this systematic 
effect is not fully understood and does not seem to be al¬ 
ways present on different lines of sight. In any case it is 
included in our estimate of the systematic error. 

Finally a small additional systematic effect is found to 
be attached to the uncertainty in the shape of the LSF, 
i.e. whether it is a single gaussian or a double gaussian 
with wide wings (indicated as “double LSF” in Table 4). 
As mentioned previously, a single gaussian LSF has a free 
FWHM in the Ht, while a double gaussian LSF has free 
amplitude ratio between both components, free FWHM 
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for each, but separation between both components fixed 
to zero pixel. The impact of the shape of the LSF on the 
measured N(D i) is shown in Figure 5 where the Ax^ is 
plotted as a function of logA^ (D i). In practice this curve 
is calculated by fixing logA^ (D i) to a given value, finding 
the best fit for this value with all other parameters free, 
and plotting the difference between this x^ and the best 
X^ obtained for all possible values of logIV(D i). In this 
figure, the curvature of a Ax^ curve gives the statistical 
error, as usual, while the relative shift between two curves 
corresponding to two different sets of assumptions is in¬ 
terpreted as an estimate of the Icr systematic error tied to 
the uncertainty on the assumptions. 

We thus conclude that the total D i column density on 
this line of sight measured using the HIRS data is: 

log A(D ^ ;l 3 39 ± 0.07 ± 0.06 

where the errors are first statistical {2a) and second sys- 
tematics (95.5% c.L). 

Combining these errors leads to: 

logIV(D i)™s ^ 13.39 ±0.09 

As can be seen from the comparison between the evalu¬ 
ations of IV (D i) made for the different apertures, there is 
no significant systematic effect tied to the aperture, con¬ 
trary to the analysis for O i and N i. The average N(D i) 
measured through the different apertures then reads: 

logIV(D i)tot = 13.40 ±0.07 

This value agrees with all previously measured values 
of IV(D i) shown in Table 4, except with that of Sahu 
et al. (1999), using the STIS-Ech#l data, which gave 
logIV(D i)= 13.55too8- appears that the dis¬ 

crepancy between this value and the others should be at¬ 
tributed to the dataset used, in agreement with the con¬ 
clusion of Vidal-Madjar (2001), and the debate around the 
value of A(D i) is now settled. 

3. THE VELOCITY STRUCTURE OF THE LINE OF SIGHT 
AND IV(H l) 

In this Section we analyze new high resolution high 
signal-to-noise STIS data of G191-B2B in order to deter¬ 
mine the number of components on the line of sight and 
provide a new estimate of the total neutral hydrogen col¬ 
umn density. As mentioned in Section 1, there exists a 
controversy in the literature about the total number of 
absorbing components, as Sahu et al. (1999) claimed to 
see two only, whereas Vidal-Madjar et al. (1998) claimed 
that three at least were present. In this Section, we present 
concrete evidence in favor of the latter. Previous measure¬ 
ments of N(H i) are compiled in Table 5; these values scat¬ 
ter around a (non-weighted) mean logV(H i)= 18.34±0.03 
(the dispersion corresponds to the non-weighted dispersion 
of the individual measurements around the mean). Indi¬ 
vidual measurements tend not to agree with each other, 
and therefore a new independent estimate is useful. More¬ 
over we argue in this Section that the total V(H i) toward 
G191-B2B is, as a matter of fact, much less well known 
than previously thought, since additional weak hot ab¬ 
sorbing components may strongly affect the column den¬ 
sity estimate from Lyman a. The new STIS data analyzed 


here represent about three times the total exposure time 
cumulated hitherto toward G191-B2B. 

All STIS observations were extracted following the 
method of Howk and Sembach (2000). We selected spec¬ 
tral lines corresponding to N i, O i. Si ii, Si ill, C ii, S ii, 
S HI in the E140H echelle configuration and Fe ii lines in 
the E230H one, both at i? ~ 90, 000 resolving power, as 
reported in Table 6. All these species are seen in one or 
more spectral lines, strong and/or weak, and cover a wide 
range of atomic masses, so that this data should provide 
strong constraints on the structure of the line of sight as 
well as on the component to component physical state and 
thermal and non-thermal broadening of each. 

In order to verify that no photospheric line is blended 
with one of the interstellar lines analyzed, we have com¬ 
puted a synthetic spectrum of G191-B2B which includes 
all species observed in the atmosphere of G191-B2B, i.e., 
G, N, O, Si, P, S, Fe, and Ni (Lanz et al. 1996; Vennes 
et al. 1996). We employ the program TLUSTY developed 
by Hubeny & Lanz (1995) to compute a NLTE metal line- 
blanketed model by using the atmospheric parameters of 
Barstow et al. (1998), i.e., logi/ = 7.4 and Teff = 54 000K, 
and the abundances determined by Barstow et al. (2001). 
From the comparison of the observations with the model 
calculation it was possible to select the lines where no or 
very weak photospheric features are coincident with the 
interstellar ones, as for the FUSE data. Furthermore it 
is possible that some photospheric lines are present but 
not predicted by the model; we thus fitted simultaneously 
as many different spectral lines of each species as avail¬ 
able, in order to identify and minimize the effect of unpre¬ 
dicted photospheric lines, as was done for the analysis of 
the FUSE data. This approach also reduces the possible 
impact of instrumental defects which may be present in 
some area of the detector but not in others. 

Following the general method of analysis described in 
the previous Section, we split the STIS data into sub¬ 
spectra, each of them centered on one spectral line and 
having typical width ~ O.SA. All sub-spectra are then 
fitted simultaneously. The wavelength zero point offsets 
are left free to vary during the fit. Note that shifts in 
the wavelength calibration from region to region could 
still be present even after a careful wavelength scale 
calibration. As an example, we detected an erroneous 
4km s“^wavelength shift between the GHRS Si ill line 
position and that in the STIS data. After comparison 
with GHRS first order grating observations, we confirmed 
that the STIS calibration was correct, and the GHRS 
calibration erroneous. We found that the average rela¬ 
tive shift between two spectral regions in the STIS data 
is ~ —O.I ± 1.1km s“^ {2a error). An rms error of 
0.5km s“^ in the STIS calibration is indeed compatible 
with the ~3.3km s“^ resolution of the instrument. The 
stellar continua in each region are interpolated by low 
order polynomials (except for Lyman a where the theo¬ 
retical continuum is used in some cases, see below), and 
the background level in each region is determined using 
the closest saturated line. The typical background flux is 
—2.6 ± 2.4% {2a error) of the continuum for the E140H 
data, and —0.10 ±4.2% {2a error) for the E230H data. Fi¬ 
nally in some cases (see below for details), the LSF is left 
free to vary and corresponds either to a simple gaussian or 
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to a double gaussian. Again, we emphasize that the uncer¬ 
tainty on these instrumental characteristics is contained in 
the final error bars given below, which we determine using 
a method, since the corresponding parameters are left 
free to vary during the fit. 

In the case of the STIS observations made with the 
E140H and E230H echelle modes, we used the tabulated 
LSE, corresponding to the slit used for these observations 
(0.2x0.2 arc sec). This is an important issue because these 
LSEs possess wings that could have some impact on the 
precise determination of the physical parameters we are 
aiming at. To test the quality of the fits we also tried 
a single gaussian LSE with free FWHM. We found ex¬ 
cellent fits with this freely variable width single gaussian 
LSE for both echelle datasets. We find for the E140H an 
average EWHM= 2.91 (±0.44) pixel (2cr) and for E230H, 
FWHM= 2.38(±0.47) pixel (2a). These determinations 
are certainly compatible with the tabulated LSE in terms 
of “average” widths, but are different in terms of shape. 
They correspond to 80,000 and 100,000 resolving power 
respectively for the E140H and E230H spectra. Finally, 
we also tested double gaussian LSEs but found no signif¬ 
icant improvement. We thus decided to present results 
with both types of LSEs (free gaussians and tabulated) in 
order to show the stability and robustness of our column 
density determinations. As we will see, such a change has 
negligible impact on the derived total H i column density. 

Following the analysis of the FUSE data, we evaluate 
the accuracy of the intrument estimate of the individual 
pixel errors by completing low order polynomial fits over 
relatively smooth and flat continuum sections close to each 
of the fitted spectral lines. These simple fits give a to¬ 
tal x^= 921.37 (sum of all iii fhe different portions of 
continuum analyzed) for 661 degrees of freedom, which 
lies well above the upper 95.5% confidence level limit of 
723.9 for a distribution with 661 d.o.f. Therefore we 
rescale all x^ by a common factor 921.37/600.63 = 1.53, 
where 600.63 corresponds to the lower 95.5% confidence 
level limit for 661 d.o.f. 

Finally, in the course of the analysis, the physical char¬ 
acteristics of each absorbing component (velocity v, tem¬ 
perature T and micro-turbulent broadening ^) are deter¬ 
mined; however we will not provide a detailed analysis of 
these characteristics and their uncertainties since we are 
primarily interested in the total number of absorbing com¬ 
ponents and the total neutral hydrogen column density. 

3.1. Number of absorbing components 

In order to settle the debate on the total number of ab¬ 
sorbing components, we have studied two and three com¬ 
ponent solutions. Simple eye inspection of the data indi¬ 
cates that at least two components are present on the line 
of sight; as mentioned previously, the red component can 
be identified with the LIC, a moderately neutral region, 
while the blue component is clearly more ionized, and, as 
we argue, more complex than a single absorber (see Vidal- 
Madjar et al. 1998 and Sahu et al. 1999). 

The result of this analysis is summarized in Table 7, 
which gives the x^ values obtained for the fits of various 
spectral regions, and for all regions (excluding Lyman a), 
for two and three component solutions, using either the 
tabulated STIS LSE or a freely varying single gaussian 


LSE. It is clear that the improvement in x^ when going 
from two to three components is extremely high whichever 
LSE is used. This is particularly true for Fe ii in which 
case the improvement can be seen directly on Figure 6. 
Indeed if one looks closely at the Fe ii lines, one can 
see an asymmetry of the line profile of the bluer com¬ 
ponent which is common to all Fe ii lines. This consti¬ 
tutes clear evidence for the presence of three absorbing 
components and not two. One sees this third component 
in Fe ii and not in other lines since Fe ii is the heaviest 
species, and its line widths are thus the smallest. More¬ 
over its lines have been observed through the echelle pre¬ 
senting the higher resolving power (100,000). This suffices 
to reveal the double structure of component B: compo¬ 
nents B1 and B2 are separated by 3.8 ± 0.5 km s“^ (2a), 
i.e. of the order of the E140H spectral resolution but 
more than that of the E230H. The radial velocities of the 
three components is estimated as ubi = 7.7 ± 0.5km s“^, 
vb 2 = 11-5 ± 0.5km s“^and wlic = 19.4 ± 0.5km s“^, with 
Icr errors due to the wavelength calibration. The LIC ve¬ 
locity is in perfect agreement with its projected velocity 
on the line of sight. 

One can also see further evidence for the double struc¬ 
ture B1-B2 in the S ill 1190.2A line shown in the upper 
right panel of Fig. 7. This line is seen only in the bluest 
component Bl, which reflects the difference in ionization 
between the various components; its spectral position is 
quite precisely known as it is located close to the Si ii line 
at 1190.4A and is seen in two different spectral orders (see 
Table 6). Note that the position of most lines is actually 
very well controled, which is important to determine the 
number of absorbing components. In effect, the sharp geo¬ 
coronal lines of N I, O I and O i* are in several locations 
(see Fig. 7), notably close to interstellar lines of N i, O i 
and Si ii. Other Si ii lines or Fe ii lines are well constrained 
since in each of these species the absorption due to the LIC 
is sharp and pronounced, and serves as a useful spectral 
reference. 

In order to test the hypothesis that the presence of a 
third component is not required by the fit, we perform an 
F—test. This test uses the Fisher-Snedecor law which de¬ 
scribes the probability distribution of the x^ ratio. Here we 
test the probability that the decrease of the x^ due to the 
inclusion of a third component could be due to the increase 
of free parameters and not to actual information contained 
in the data. The result is shown species by species and for 
all species fitted simultaneously in Table 7, and gives a 
probability < 10“^ that this third component is not re¬ 
quired by the data. The F—tests performed species by 
species also confirm our previous impression that Si ii and 
Fe II are the two most sensitive species to the presence of 
this component. We also note that the need for a third 
component is present whichever LSF is used. 

We note that the x^ values corresponding to the use of 
the tabulated STIS LSF are too high, as they all lie above 
the upper 95.5% confidence level limit for statistical fluc¬ 
tuation (an except being Fe ii for three components). In 
contrast, the x^ values for a freely varying single gaussian 
LSF are satisfactory, since only O i and Si ii stand out with 
X^ beyond the 95.5% confidence level for the three com¬ 
ponent solution. Their high x^ propagates into the total 
X^ summed over all windows. However the Si ii 1190A, 
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1193A and O i 1302A lines are saturated, and small errors 
in the estimates of the error bars at the bottom of the 
line could potentially explain large differences in x^- 
effect our procedure of rescaling the amounts to rescal¬ 
ing all error bars by a common factor, which may not 
be correct for pixels at the bottom of saturated lines for 
the following reason. Error bars on flux values arise as a 
combination of background noise, fixed pattern noise and 
Poisson noise. Background noise dominates for pixels with 
flux values close to zero, while for high S/N data Poisson 
noise or fixed pattern noise would dominate in the contin¬ 
uum. If the errors in the estimation of the noise array for 
the background component and for the continuum com¬ 
ponent (fixed pattern or Poisson) are different, then the 
procedure of rescaling error bars with a common factor is, 
strictly speaking, incorrect. Finally, since this rescaling 
factor was estimated from the continuum, this rescaling 
should be correct for all pixels whose error is not dom¬ 
inated by the other noise component (background). We 
also note that this O i line has a complex profile, with two 
geocoronal and one nearby photospheric feature (see Fig¬ 
ure 7). Nevertheless, we decided to keep this line in our 
analysis as it adds important constraints on the B1 and 
Lie component T and ^ parameter evaluations. 

Using the F—test, we can determine the probability that 
the use of the freely varying Gaussian LSF is not required 
(instead of using the STIS LSF). In all cases, it is found 
that this probability is smaller than 10“"^, and therefore we 
conclude that the gaussian LSF with a width that varies 
from spectral region to spectral region is required to fit 
the STIS data. This result may be unexpected but it can 
find a simple explanation in the fact that the tabulated 
STIS LSF is an average LSF given over the whole spectral 
range, while our gaussian LSF has a FWHM that varies 
with wavelength. For standard datasets it is probably suf¬ 
ficient to use the STIS LSF but the very high quality of the 
present data necessitates higher order corrections to this 
LSF. In particular we noted a slight broadening of the LSF 
near the order edges. Furthermore, the wings of the tab¬ 
ulated LSF are not so important and the single gaussian 
LSFs are sufficient to properly fit the data. 

We thus conclude that three absorbers are present on 
the line of sight to G191-B2B. 

3.2. Lyman a 

The background level is directly measured at the bot¬ 
tom of the strongly saturated Lyman a line. It is found 
to be of the order of —1.0(±1.0)% {2a) of the nearby con¬ 
tinuum, i.e. within a few percent. Variations in the pho¬ 
tospheric continuum over the interstellar Lyman a line 
presents an additional difficulty. Indeed the instrument 
sensitivity over this spectral region is evaluated by using 
model spectra of white dwarfs and comparing them to the 
data, and one of the calibration WDs used is G191-B2B 
itself. Therefore it is difficult to separate the instrumen¬ 
tal effects from the photospheric intrinsic profile in the 
spectrum processed through the STIS pipeline. For this 
reason, we follow the approach detailed in Vidal-Madjar et 
al. (1998) and further discussed by Vidal-Madjar (2000), 
and evaluate simultaneously all parameters and the stel¬ 
lar continuum, which we model in two different ways. In 
the first approach, we model this continuum as a poly¬ 


nomial whose parameters are adjusted during the fitting 
procedure, and which represents the real photospheric pro¬ 
file times the instrument sensitivity; we denote this fitting 
model “ L”, which stands for uncorrected (as the data have 
not been corrected prior to the fit by a stellar model). In 
the second approach, we correct the data prior to the fit 
by a stellar model and keep a free low-order polynomial 
during the fit which models the inaccuracy of the stellar 
model and the variation of the instrument sensitivity; we 
denote this approach by “ C”. The calculated stellar model 
was shifted by 24.56 km s“^ in order to be consistent with 
the velocity of the photospheric features clearly detected 
in N V and Si ill (see Figure 7). These two approaches have 
been followed in parallel to estimate the uncertainty in the 
modeling of the stellar continuum on the final solution. 

We use two models to compute the photospheric Ly¬ 
man a line profile. First, we use the program SYNSPEG 
(1. Hubeny 2000, private communication) in conjunction 
with the TLUSTY NLTE metal line-blanketed model de¬ 
scribed above. This version of SYNSPEG contains the 
Lemke’s Stark broadening tables for hydrogen, which were 
computed by Lemke (1997) within the framework of Vi¬ 
dal, Gooper, & Smith (1973). Second, we use Detlev 
Koester’s LTE code and his best fit parameters derived 
under the assumption of pure-H LTE (Tea = 60880 K and 
log(/ = 7.59). Figure 8 shows the NLTE and LTE Lyman 
a line profiles and illustrates that NLTE effects are signif¬ 
icant only in the core of the line, which is formed high in 
the atmosphere where departures from LTE are important 
(Wesemael et al. 1980). However, as is shown in Fig. 8 
the observations do not contain any information in that 
central region of the photospheric line simply because it is 
lost at the bottom of the saturated interstellar Lyman a 
line. Our fitting procedure will however be able to test the 
slightly different slopes in the wings of the models. 

An additional correction made by a low-order polyno¬ 
mial is needed not only to take into account uncertainties 
in the instrument sensitivity but also because when one 
fits a section of the G191-B2B stellar continuum devoid of 
stellar or interstellar absorption, it is necessary to model 
the continuum with a polynomial of degree of order 1 to 3 
for a spectral width of order 0.3A (see for instance Fig. 7). 
For a spectral region of width ~ 3A as is the case for 
Lyman a, a 6*^ order polynomial provides a satisfactory 
approximation to the stellar continuum. 

The best-fit solutions obtained for 3 components on the 
line of sight for all lines in the STIS domain including Ly¬ 
man Of and for various models (prior correction or not of 
the Lyman a stellar continuum, type of LSF) are summa¬ 
rized in Table 8. The zero point wavelength of the Lyman 
a spectral domain was found to be compatible with the 
other regions well within the Icr uncertainty of 0.5km s“^. 

Goncerning the fit with a free double gaussian (model 
5 in Table 8), we found that the wings of the broader 
gaussian (about 5 to 10 times larger than the narrow one) 
contribute in amplitude (relatively to the amplitude of the 
narrow core gaussian) to less than 1% in all but 5 spec¬ 
tral windows, less than 2% in 3 of these 5 and to about 
3% in the 2 remaining ones. Since this effect is marginal, 
we kept only as comparison tests the use of both the tab¬ 
ulated STIS LSF or the simple gaussian LSF. All these 
effects were included in the evaluation of the statistical 



errors. 

It appears that different values of A^(H i) are obtained 
depending on the model used; in particular the comparison 
of models 1 and 2 on the one hand, and of models 3 and 4 
on the other hand shows that the value of A^(H i) derived 
depends whether the stellar continuum at Lyman a has 
been corrected or not beforehand by a theoretical stellar 
profile. One obtains logiV(H i)=18.37 for an uncorrected 
profile, and logiV(H i)=18.32 for a corrected stellar profile. 
Since there is no significant difference in between these 
solutions, this discrepancy must be interpreted as a sys¬ 
tematic uncertainty tied to the choice of modeling of the 
continuum. A Ay^ analysis around each of these solutions 
give a statistical error of ±0.01 dex (Icr), so that the value 
of N (H i) is dominated by the above systematic error. Fur¬ 
thermore after correction of the stellar continuum by an 
LTE calculation, which as shown in Fig. 8 represents an 
intermediate situation between NLTE and non-corrected 
stellar profiles, the value of N(H i) derived was found to 
be an intermediate value between the above two. 

However a further systematic error appears if additional 
weak hot components are present on the line of sight, and 
its investigation is the subject of the following Section. 

3.3. Systematic uncertainties and additional hot 
components 

One cannot exclude the presence of weak H i compo¬ 
nents which could perturb or bias the measurement of 
N(R i), but could not be detected in any other species 
due to their weak column density. Such absorbers could 
arise as the signature of high velocity shocks (Cowie et al. 
1979), or as cloud interfaces with the hot gas within the 
local ISM (Bertin et al. 1995) or as “hydrogen walls”, i.e., 
the shock interaction between the solar wind (or stellar 
wind) and the surrounding ISM (Linsky 1998). The lat¬ 
ter possibility has been modeled by Wood et al. (2000), 
and a prediction for the line of sight toward G191-B2B is 
shown in Fig. 9 overlaid on the STIS data. It shows that 
most of the expected absorption should take place in the 
saturated core of the interstellar line, but some weak ab¬ 
sorption (~ 5%) may be present, extending over several 
tenths of an Angstrom on the red side of the line, due to 
the neutral hydrogen atoms seen behind the shock in the 
downwind direction where G191-B2B is located. 

We have thus investigated the possible impact of ad¬ 
ditional weak hot components on the determination of 
N(H i). In order to do so, we have added one or two addi¬ 
tional components in H i only, and performed the fit of all 
lines as before with an NLTE correction of the Lyman a 
profile and using a free single Gaussian LSF (model 3). In 
order to constrain the presence of these extra absorbers, 
we have added to the fit the four Lyman (3 lines observed 
through the FUSE HIRS aperture. For each new model 
characterized by the number of additional components, we 
have measured A(H i) using a Ay^ analysis. The best fit 
solutions with zero, one and two extra absorbers are shown 
for Lyman a and Lyman /3 in Fig. 10. We have found 
that this introduction of additional components leads to 
a substantial decrease of the total A(H i) together with 
a significant improvement of y^. In particular, we found 
that the best fit solution for one extra component leads 
to Ay^ = 26 (1961 d.o.f.) and logA(H i) = 18.24, and 


two extra components lead to Ay^ = 39.4 (1958 d.o.f.) 
and logA(H i)=18.11. The y^ values quoted have been 
rescaled by a factor 1.53 corresponding to the STIS data 
obtained in the previous section, and which remains close 
(20%) to the factor measured in the FUSE data range. 
Moreover, if a Ay^ analysis is performed for each model, 
the curve obtained flattens for low values of N (H i) down 
to logA(H i)=18.0; in other words, the space of solutions 
becomes degenerate in the low A(H i) region, and con¬ 
sequently the value of A(H i) measured from the profile 
fitting of the Lyman a region is subject to a large sys¬ 
tematic uncertainty. This behaviour can be explained as 
follows. In the three component best fit solution, most 
of the hydrogen is contained in the LIC whose absorption 
makes the red wing of Lyman a. In the solutions with 
extra absorbers, the combined profile of these hot compo¬ 
nents can replace the damping wing contribution of the 
LIG, as shown in Fig. 10. 

In order to scan the parameter space, the A(H i)lic has 
been forced to take various values, since the profile fitting 
code cannot force the total H i column density to take a 
particular value. When the LIC A(H i) becomes negligi¬ 
ble with regards to the total A(H i), we force the next 
dominant component to take various values of A(H i). In 
this way it is possible to scan the parameter space of the 
total A(H i). One finds that once the LIC has become 
negligible, component B2 makes most of the total A(H i), 
and that when the total A(H i) is taken below 18.0, the 
y^ rapidly increases, which gives the final error bar on the 
total A(H i) estimate. 

At this stage, one should point out that these best-fit so¬ 
lutions are quite difficult to find. The impact of additional 
weak absorbers on the profile fitting of N (H i) had already 
been studied by Vidal-Madjar et al. (1998) but no signif¬ 
icant effect had been found, as the above solutions had 
eluded detection. The above solutions have actually been 
obtained only in a late stage of the present work. We in¬ 
terpret this as evidence for the fact that the y^ surface be¬ 
comes complex when additional unconstrained absorbers 
are introduced in the fit of the fully blended Lyman a 
profile and optimization is then delicate. 

One should also remark that we have modeled these 
extra absorbers using Voigt profiles, which implicitly as¬ 
sumes a Maxwellian distribution for the velocities of the 
atoms. However if these absorbers correspond to interstel¬ 
lar structures similar to hydrogen walls, the latter assump¬ 
tion is incorrect and the overall profile should be closer to 
that modeled by Wood et al. (2000) and shown in Fig. 9. 
For this reason, the above modeling and solutions should 
be interpreted as a first approximation. For this reason 
as well, we do not discard the values for A(H i) obtained 
in the previous section without additional absorbers, even 
though the presence of extra absorbers induce a significant 
gain in y^. At the worst the final value for A(H i) will be 
more conservative. With these remarks in mind, we con¬ 
clude that the total hydrogen content toward G191-B2B 
should be contained in the interval: 

18.00 < logA(H i) < 18.37. 

In the absence of any indication on the distribution of 
the errors for A(H i), we consider the above two extreme 
values as 2a limits since they correspond to rather ex- 
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treme solutions, and quote the following value for A^(H i): 
log-/V(H i)= 18.18 ± 0.09 (Icr error). One should note 
that this effect of additional absorbers had not been found 
in previous studies of the total hydrogen content toward 
G191-B2B. Therefore it is important to remark that the 
above large uncertainty must also affect the previous val¬ 
ues of A^(H i). In other words, the above results supersedes 
any previous estimate of 7V(H i) made toward G191-B2B 
using profile fitting of the Lyman a line. One may also 
wonder whether extra absorbers may affect the determi¬ 
nation of fV(H i) toward other stars; this important ques¬ 
tion is the subject of a forthcoming paper (Vidal-Madjar 
& Ferlet 2002). On the other hand, the above value for 
fV(H i) should be contrasted with the value derived from 
EUVE observations and modeling of the atmosphere of 
G191-B2B which gave log7V(H i)=18.315± 0.013 (Dupuis 
et al. 1995, 2a error), log7V(H i)=18.32 (Lanz et al. 1996, 
no error bar quoted) and the recent detailed measurement 
logA^(H i)=18.30 ± 0.09 (Barstow & Hubeny 1998; this 
value also agrees with the more recent work of Barstow, 
Hubeny & Holberg 1999). These values would tend to 
indicate that the contamination of the Lyman a profile 
by weak hot absorbers (if any) is not important, as they 
agree with the values we obtained without including such 
extra absorbers. However, as noted in Barstow & Hubeny 
(1998), these fits of the EUVE spectrum of G191-B2B typ¬ 
ically lead to very large and unexplained reduced im¬ 
plying that the overall fit is not yet satisfactory in spite of 
the use of sophisticated WD atmosphere models, and that 
some unknown effects have yet to be accounted for. 

4. DISCUSSION AND CONCLUSIONS 

We have measured the total column densities of D i, N i 
and O I toward G191-B2B using unsaturated absorption 
lines of these elements in high quality FUSE spectra. After 
a careful examination of the possible systematic uncertain¬ 
ties tied to the choice of the stellar continuum and to the 
instrumental configuration, we have derived the following 
column densities with 2a uncertainties: 

logV(D i)tot = 13.40 ±0.07, 
logiV(0 i)tot = 14.86 ±0.07, 
logV(N i)tot = 13.87 ±0.07. 

We have also analyzed new high signal-to-noise ratio 
high resolution STIS observations of G191-B2B and pro¬ 
vided concrete evidence for the presence of at least three 
interstellar absorbing components on the line of sight by 
analyzing the interstellar absorption lines of N i, O i. Si ii. 
Si III, S III and Fe ii present in the STIS bandpass. We 
have also measured the total hydrogen column density on 
the line of sight using the velocity structure derived from 
the above metals. We have performed an exhaustive study 
of systematic effects on the value of V(H i). In particular 
we have discovered a new major source of uncertainty on 
V(H i) tied to the possible presence of additional weak 
hot absorbers whose combined absorption profile can con¬ 
tribute significantly to the wings of the blended Lyman a. 
profile. The column density of these absorbers is small 


compared to the other main components, and they would 
not be detected in any other species than H i, but their 
contribution to the Lyman a. absorption profile can re¬ 
duce significantly the total H i column density measured 
from the profile fitting. In order to constrain their impact, 
we have analyzed simultaneously Lyman a and the higher 
order Lyman lines, and concluded that the best value of 
V(H i) toward G191-B2B is: 

log V(H i)tot = 18.18 ± 0.18 (2(jerror) 

We emphasize that this uncertainty is a systematic un¬ 
certainty which had gone unnoticed before. Therefore the 
above result supersedes previous estimates of V(H i) to¬ 
ward G191-B2B obtained from the profile fitting of Ly¬ 
man a. A detailed analysis of this uncertainty and its 
consequences on V(H i) determinations toward other stars 
is discussed in a companion paper (Vidal-Madjar & Fer¬ 
let 2002). We thus derive the following neutral abundance 
ratios toward G191-B2B, with 2a uncertainties: 


(D/H)tot = 1 . 66 l°;® 10 -® 
(D/O)tot = 3.49 ±0.7810"^ 
(D/N)tot = 3.41 ±0.7610-1 
(0/H)tot = 4.79l^;^ lO'^ 
(N/H)tot = 4.9012-6 10-5 


Most of the uncertainty in the above result results from 
the systematic uncertainty on the V(H i) determination. 
This clearly shows the importance of measuring accurate 
(D/0) and (D/N) ratios in the interstellar medium in¬ 
stead of abundances relative to hydrogen, as emphasized 
by Timmes et al. (1997). Interestingly if one uses the 
recent measurement of V(H i) from the modeling of the 
atmosphere of G191-B2B and the fit of the EUVE spec¬ 
trum, log V(H i)= 18.30 ±0.09 {2a), one finds (D/H)tot = 
1 .261 q;29 ^ 40“® (2fr). At this stage, however, due to the 
uncertainty inherent to the modeling of the white dwarf 
atmosphere, it is probably more conservative to use the 
interstellar determination for V(H i), and therefore the 
previous value of the (D/H) ratio. 

The above new value for the (D/H) ratio agree with the 
range of values measured by Linsky (1998) toward a dozen 
stars of the LISM and with the values previously derived 
toward G191-B2B. However the discrepancy between pre¬ 
vious estimates of the (D/H) ratio toward G191-B2B and 
the LISM average D/H ratio has disappeared due to a re¬ 
vision of the uncertainty on the estimation of the total H i 
content. A detailed interpretation of this (D/H) value and 
of the accompanying (D/0) and (D/N) ratios and their 
implications is provided in a companion paper by Moos et 
al. ( 2002 ). 
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Table 1 


Strongest N i and O i spectral lines in the FUSE domain. 


Wavelength (A) 

Element 

/ 

^ (s-i) 

Comment 

1025.4440 

D I 

0.264 10-1 

0.190 10^ 


1025.4429 

D I 

0.527 10-1 

0.190 10® 


972.2725 

D I 

0.967 10-2 

0.813 10® 


972.2721 

D I 

0.193 10-1 

0.813 10® 


949.4848 

D I 

0.465 10-2 

0.421 10® 

Weak 

949.4846 

D I 

0.929 10-2 

0.421 10® 

Weak 

1134.9803 

N I 

0.435 10-1 

0.150 10® 


1134.4149 

N I 

0.297 10-1 

0.154 10® 


1134.1653 

N I 

0.152 10-1 

0.158 10® 


954.1042 

N I 

0.676 10-2 

0.330 10® 

Blend 

953.9699 

N I 

0.348 10-1 

0.170 10® 


953.6549 

N I 

0.250 10-1 

0.183 10® 


953.4152 

N I 

0.132 10-1 

0.193 10® 


1039.2303 

0 I 

0.920 10-2 

0.947 10® 


1026.4757 

0 I 

0.246 10-2 

0.111 10 ® 

Weak 

1026.4744 

0 I 

0.187 10-3 

0.118 10 ^ 

Weak 

1025.7633 

0 I 

0.201 10-3 

0.212 10 ^ 

Blend 

1025.7626 

0 I 

0.302 10-2 

0.191 10® 

Blend 

1025.7616 

0 I 

0.169 10-1 

0.765 10® 

Blend 

988.7734 

0 I 

0.465 10-1 

0.226 10 ® 

Strong 

988.6549 

0 I 

0.830 10-2 

0.566 10® 

Blend 

988.5778 

0 I 

0.553 10-3 

0.629 10^ 

Weak 

976.4481 

0 I 

0.331 10-2 

0.386 10® 


971.7382 

0 I 

0.116 10-1 

0.585 10® 

Strong 

971.7376 

0 I 

0.207 10-2 

0.146 10® 

Blend 

950.8846 

0 I 

0.158 10-2 

0.194 10® 


948.6855 

0 I 

0.631 10-2 

0.100 10 ® 


936.6295 

0 I 

0.365 10-2 

0.100 10 ® 


929.5168 

0 I 

0.229 10-2 

0.100 10 ® 


925.4460 

0 I 

0.354 10-3 

0.459 10^ 

Weak 

924.9500 

0 I 

0.154 10-2 

0.100 10 ® 


921.8570 

0 I 

0.100 10-2 

0.562 10^ 

Weak 


Table 2 

FUSE OBSERVATION LOG OF G191-B2B. 


Dataset 

Aperture 

^exp (kS6c) 

iVexp 

Date 

S3070101 

LWRS 

15.5 

32 

2000.01.14 

P1041202 

MDRS 

15.5 

21 

2000.01.13 

P1041201 

HIRS 

15.5 

32 

2000.11.06 




10”' 'erg/cm^s/A 
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920 940 960 980 

Wavelength (A) 


Fig. 1.— FUSE HIRS spectrum of G191-B2B (segment SiClB). The Lyman series is clearly seen down to the Lyman limit. The positions 
of the D I, O I and N I lines considered in our study are indicated. 


Table 3 

Published and FUSE O i and N i column densities. 


Spectrograph / 
Aperture 

logiV(0 i)tot 
2(7 error 

log07(N i)tot 
2(7 error 

Reference 

HST GHRS-Ech. 

14.84±0.04 

13.90±0.02 

Vidal-Madjar et al. 1998'^ 

FUSE LWRS+MDRS 

14.79±0.04 

13.82±0.07 

This work‘d 

FUSE MDRS 

14.84±0.08 

13.89±0.06 

This work^ 

FUSE HIRS 

14.88±0.06 

13.84±0.07 

This work*’ 

EUSEUkW 

14.86±0.07 

13.87±0.07 

This work*’ 


^'Three components were considered in the fit 
'^One interstellar component was considered in the fit 
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971.5 972 972.5 973 


Wavelength (A) 


Fig. 2.— FUSE MDRS spectrum of G191-B2B around Lyman 7 (segment SiC2A). The solid line in each panel shows the best fit obtained 
by a simultaneous profile fitting of all H i, D i, O i and N i lines of all FUSE datasets. Clearly these observations do not resolve the multiple 
absorbers along the line of sight. In the upper panel the continuum is interpolated by a smooth polynomial, and three interstellar components 
corresponding to the velocity structure of the line of sight derived from the higher resolution STIS observations (see Section 3) were considered; 
in the middle panel, the stellar continuum is also interpolated by a polynomial, while only one absorbing component is considered. In the 
lower panel the stellar continuum has been corrected by a theoretical NLTE stellar profile, the residual continuum is modeled by a polynomial 
and only one absorbing component has been assumed. See text for details. 
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953 953.5 954 954.5 


Wavelength (A) 


Fig. 3. — FUSE spectra of G191-B2B around the N i multiplet at 954A. The aperture, segment and resolving power R corresponding to each 
dataset are indicated on the figure. The line at 954.1 A was excluded from the profile fitting since it is likely contaminated by a photospheric 
feature, as indicated by theoretical modeling of the stellar continuum, and as shown clearly in the upper and middle panels; in the lower 
panel, the apparent absence of this line may result from a statistical fluctuation or from a detector feature. The higher resolution of the HIRS 
aperture data is apparent. 
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Fig. 4.— FUSE HIRS spectra of Lyman (3 and Lyman 7 , recorded on the segments indicated in the figure. The solid line shows the final 
best-fit solution to all lines of H i, D 1 , N 1 and O 1 in all HIRS data fitted simultaneously. The weak emission line at the bottom of the Lyman 
lines is due to H 1 geocoronal emission; the pixels affected by this emission are not considered in the fit. The stellar continuum is interpolated 
by a smooth polynomial. See text for details. 
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Fig. 5. — Curve of deviations around the best fit a function of logA^(D i). Only the curve for the HIRS data is shown here; 

statistical error bars for other aperture data were obtained using the same method. The Ax^ values shown here have been rescaled to 
compensate for uncertainties in the individual pixel errors (see text). The various curves correspond to different models for the fit: in thick 
line, the D i and H i are fitted simultaneously with a LSF modeled as a single gaussian with free FWHM; in intermediate thickness, the same 
approach but with a double gaussian LSF, with free amplitude ratio and FWHMs; in thin line, H i is excluded from the fit, the continuum to 
the D I absorption is modeled by a polynomial, and the LSF is a simple gaussian. The curvature of these curves give the the statistical errors 
while their relative shifts give an estimate of the overall systematic uncertainty related to the different models. The hatched area shows the 
final 95.5% confidence level error for logA^(D i) measured using the HIRS data. 
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Fig. 6 .— The STIS data covering the three Fe ii lines observed with the E230H Echelle grating. All fits shown here were made using a 
freely varying single gaussian LSF; the Fe ii lines were fitted simultaneously with all other species including Lyman a (see text). The left 
panels show the fits with two absorbing components only, and the right panels show the fits using three components. One can clearly see the 
same asymmetry of the bluer component B in all Fe ii lines which reveals the presence of complex substructure in this component, hence the 
need for 3 absorbers in total. The corresponding to these fits is 227.9/139 for the two component solution, and 123.3/136 for the three 
component solution; the gain is clear and comforts the visual impression. The error bars on each pixel are of order of 2% of the continuum 
level. 






18 



-50 0 50 -50 0 50 


V (km/s) V (km/s) 


Fig. 7.— The STIS data with the best-fit three component solution using the tabulated STIS LSF; only 8 sub-spectra out of 19 in total 
are shown here, and are labeled with the central spectral line wavelength and Echelle order: the N i 1200A triplet along with the nearby 
geocoronal absorptions (marked “Geo.”, in two different locations because two data sets taken at two different epochs are here averaged), the 
O I line with the corresponding geocoronal O i absorptions and a nearby photospheric feature (noted “Phot.”), three of the Si ii lines with 
either the nearby S ill feature or the O i* geocoronal absorption and the Si ill line with an additional broad Si ill photospheric line located at 
the same velocity shift as the N v photospheric line (which was fitted simultaneously). Note also that the S ill lines clearly detected in two 
different orders near 1190.2A are also spectrally well located. The fact that they show up only in the bluest component B1 is an additional 
reason why 3 and not 2 components are needed along that line of sight. 
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Fig. 8. — The STIS data over the Lyman a spectral region along with the LTE (dashed line) and NLTE (dotted lines) photospheric profiles. 
Three different NLTE profiles for atmospheric parameters in the ranges 54000 K < Tgff < 55000 K and 7.5 < \ogg < 7.6 are shown. The 
LTE calculation uses Teff = 60880 K and log^ = 7.59. The radial velocity for the photospheric profiles is set to +24.56km s“^. The two thin 
vertical lines delimit the wavelength range in which no information on the stellar profile is contained in the data. The fitting procedure thus 
tests only the difference between the wings of the models and not their cores whose difference is more pronounced. 



1214.5 1215.0 1215.5 1216.0 1216.5 1217.0 

Wavelength [A] 


Fig. 9.— The Lyman a interstellar absorbtion observed with STIS with the hydrogen wall model of Wood et al. (2000) for the line of sight 
of G191-B2B overplotted in dashed line. The STIS data has been normalized by the best-fit continuum. 
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Fig. 10.— Best-fit solutions for all lines in the STIS data plus four Lyman (3 from the FUSE HIRS dataset using three main interstellar 
components plus zero, one and two extra H i components (upper, middle, and lower panels respectively). The total absorption profile is shown 
as the solid line, and the dotted and dashed lines show the individual profiles (convolved with the LSF); the thick dashed lines in particular 
show the contribution of the extra absorbers. The solutions are labeled with the total rescaled (foi^ 1964, 1961, and 1958 d.o.f. respectively 
for the upper, middle and lower panels) and total H i column density (in log). 
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Table 4 

Published and FUSE D i column densities. 


Spectrograph / 
Spect./Aper. 

log A^(D i)tot 
2(7 error 

Reference 

/Model 

HST GHRS-Ech. 

13.43±0.02 

Vidal-Madjar et al. 1998 ^ 

HST STIS-Ech.#l 

10.00_Q_Q^ 

Sahu et al. 1999 

HST GHRS-Ech. 

13.40±0.04 

Sahu et al. 1999 

FUSE LWRS-fMDRS 


D I, no H I, this work 

FUSE MDRS 

13.41±0.09 

D I & H I, this work ^ 

EUSE HIRS 

13.36±0.08 

D I & H I (Fig. 4), this work 

FUSE HIRS 

13.38±0.06 

D I & H I, double LSF, this work ® 

FUSE HIRS 

13.42±0.08 

D I, no H I, this work ® 

EUSE All 

13.40±0.07 

this work * 


free absorbing components assumed in the profile fitting, and stellar con¬ 
tinuum modeled with a low order polynomial free to vary during the fit 

^2 free absorbing components in the profile fitting, and stellar continuum 
fixed and modeled by NLTE calculations 

'^1 free absorbing component, stellar continuum modeled by a freely varying 
low order polynomial, and D i fitted alone without the H i line 

'^1 free absorbing component, stellar continuum modeled by a freely varying 
low order polynomial, and D i fitted simultaneously with the H i line 

®1 free absorbing component, stellar continuum modeled by a freely varying 
low order polynomial, and D i fitted simultaneously with the H i line, LSF 
profile modeled with a double Gaussian 

^combination of FUSE models c, d and e 


Table 5 

Published H i column densities. 


log N{H i)tot (2cr) 

Spectrograph 

# Comp. 

Continuum 

Reference 

18.315±0.013 

EUVE 

— 

— 

Dupuis et al. 1995 

18.32 - 

EUVE 

— 

— 

Lanz et al. 1996 

18.36±0.04 

GHRS-G160M 

3 

free 

Lemoine et al. 1996 

18.38±0.02 

GHRS-Ech. 

3 

free 

Vidal-Madjar et al. 1998 

18.30±0.09 

EUVE 

— 

— 

Barstow & Hubeny 1998 

18.31±0.03 

STIS-Ech. #1 

2 

fixed 

Sahu et al. 1999 

18.34±0.02 

GHRS-Ech. 

2 

fixed 

Sahu et al. 1999 



Table 6 

Spectral lines in the STIS domain used for the line of sight structure study. 


Wavelength (A) 

Species 

/ 

Al (s-l) 

Ech. Order 

1334.5320 

C II 

0.128 

0.288 10® 

316 

1200.7098 

N I 

0.0430 

0.398 10® 

350, 351 

1200.2233 

N I 

0.0862 

0.399 10® 

350, 351 

1199.5496 

N I 

0.130 

0.401 10® 

351 

1242.8040 

N V 

0.500 

0.336 10® 

339 

1302.1685 

0 I 

0.0519 

0.340 10® 

323 

1304.3702 

Si II 

0.0917 

0.107 10^® 

323 

1193.2897 

Si II 

0.585 

0.409 lO^o 

352, 353 

1190.4158 

Si II 

0.293 

0.410 10^0 

353, 354 

1206.5000 

Si III 

1.67 

0.255 10^® 

349 

1259.5190 

S II 

0.0166 

0.465 10“ 

334 

1190.2030 

S III 

0.0231 

0.651 10® 

353, 354 

2382.7651 

Fe II 

0.320 

0.313 10® 

324 

2586.6499 

Fe II 

0.0691 

0.272 10® 

298, 299 

2600.1729 

Fe II 

0.239 

0.270 10® 

297 


Table 7 

^ comparison for the two and three component solutions using either the tabulated STIS LSF OR freely 

VARYING SINGLE GAUSSIAN LSF. 


Species 

X® / d.o.f. 

2 comp. 
STIS LSF 

^ ! d.o.f. 

3 comp. 
STIS LSF 

F—test prob. 

2 Rs 3 comp. 

/ d.o.f. 

2 comp, 
free Gaussian 

X^ / d.o.f. 

3 comp, 
free Gaussian 

F—test prob. 

2 vs i comp. 

N I 

395.0/314 

382.3/311 

1.7% 

323.2/309 

309.1/306 

0.3% 

0 I 

214.3/61 

193.1/58 

11% 

113.9/60 

101.4/57 

8.3% 

Si II 

772.0/468 

651.4/465 

< 0.01% 

678.6/463 

553.6/460 

< 0.01% 

Si iiF 

210.6/152 

192.8/149 

0.4% 

185.2/150 

170.8/147 

0.75% 

Fe II 

277.6/142 

141.8/139 

< 0.01% 

227.9/139 

123.3/136 

< 0.01% 


2677.8/1270 

1908.8/1259 

< 0.01% 

2076.7/1250 

1442.5/1239 

< 0.01% 


^^the has been rescaled to compensate for the inaccuracy of the noise array (see text) 

*’The F—test gives the probability that a third absorbing component is not required by the model 
'^Fitted simultaneously with the N v line to control the photospheric Si ill line 
‘^Includes C ii, S ii and S ill 




Table 8 


Three component best-fit solutions for all lines in the STIS region including Lyman 


Fit 

model 

VBl 

(km s“^) 

fB2 

(km s“^) 

^'LIC 

(km s“^) 

logiV(H i)tot 

1753 d.o.f. 

^a,c 

7.7 

11.2 

19.4 

18.32 

2469.9 

2b,c 

7.7 

11.3 

19.4 

18.37 

2470.1 

3^3- 

7.7 

11.7 

19.4 

18.32 

2006.8 

4b,d 

7.7 

11.6 

19.4 

18.37 

2001.8 

ga.e 

7.7 

11.6 

19.4 

18.33 

1987.3 


'^data normalized beforehand by NLTE stellar continuum, and stellar con¬ 
tinuum residuals modeled during the fit by a 6**' order polynomial 

unnormalized data, with stellar continuum modeled during the fit by a 6*^ 
order polynomial 

‘^LSF corresponds to the tabulated STIS LSF 

'^LSF profile modeled with a single Gaussian (free to vary during the fit) 
®LSF profile modeled with a double Gaussian (free to vary during the fit) 



